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\ Abstract. The abihty to measure characteristics of source shapes using non- 

' identical particle correlations is discussed. Both strong-interaction induced and 

Coulomb induced correlations are shown to provide sensitivity to source shapes. 
By decomposing correlation functions with spherical or Cartesian harmonics, 
\ details of the shapes can be especially well isolated. 
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|n ' The space-time development of the expanding fireball in a heavy-ion collision 

, at RHIC is driven by the pressure and viscosity of the matter during the novel 

' stage of the reaction. Insight into these fundamental questions can be gained by 

understanding the spatio-temporal aspects of the collision. Correlation functions [ 
Q ' n El El provide a direct link to the space-time structure of the source function, 

^\ 5(P,r), 



i?(P, q) = 1 d\' [|0(q', r')|' - l] 5(P, r'). (1) 



' 

^ , Here, i?(P, q) + 1 is the ratio of two-particle cross-sections to the corresponding 

\^ ' quantity for mixed-events, and is zero for uncorrelated emission. The source func- 

5^ , tion ^(P, r') measures the probability that two particles of the same velocity whose 

total momentum is P are separated by r. The primes denote that the coordinates 
are measured in the rest frame of the pair. By measuring i?(P,q) as a function 
of the relative momentum q' = (pjj — Ph)/2, one can exploit the structure of the 
relative wave function 0(q',r') to extract information about S'(P,r). 

Understanding the source function in Eq. assists in the determination of 
collective flow and lifetimes which directly address questions concerning the equa- 
tion of state. In this talk the phenomenology that connects the source function 
to these important issues will be ignored. Instead, this talk focuses on the ability 
to determine size and shape information about S'(r') from i?(q'). Thus, from this 
point on the total momentum label and the primes are suppressed. This problem 
can be considered as an inversion problem. Both i?(q) and S{y) can be considered 
as vectors and |</)(q, r)p can be considered as a matrix. Simply stated. Given i?(q) 
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Fig. 1. Coulomb repulsion suppresses trajectories where the outgoing relative 
momenta is anti-parallel to the original separation. For sources biased to one side, 
as illustrated in the upper panel, the correlation function is suppressed for opposite 
directions. For sources with a quadrapole distortion along a given axis, the distor- 
tion is strongest for emission parallel or anti-parallel to the axis as illustrated in the 
lower panel. 



and |0(q, r)p, what can one determine about S'(r)? For non-interacting identical 
particles, |0p = l + cos2q-r, and one can Fourier transform R{q) to find S{r). Our 
goal is to investigate the analogous deconvolution for wave functions whose forms 
are driven by Coulomb and strong interactions. 

A straight-forward case to consider is interactions through classical Coulomb 
correlations. In the classical limit, 

\(j){q, r, cos 6',r)p 



where qo is the relative momentum when the particles were separated by r, 9qr is 
the angle between q and r, and x — 2/ie^/g^r. When integrating over angles, |</>p = 
Vl — X, and the correlation function approaches unity at large g as 1 — {1 / r) ^e^ / . 
By working inward from large the correlation function can be used to determine 
inverse moments of the source function, 1/r, etc. 

The cos 9qr dependence in Eq. (jSJ can be exploited to extract shape informa- 
tion about the source. This sensitivity is driven by the scattering of the Coulomb 
trajectories which lowers the population for angles where q leaves anti-parallel to r 
as illustrated in Fig. 

The quantum Coulomb wave function is more complex, and adds a additional 
sensitivity to qr/Ti in addition to x and cos 9 qr- Quantum effects tend to smear out 



d?q 



(2) 



1 -I- cos 9, 



qr 



■v/(l -I- COS 9qr — x)^ — X- 



:e(l -I- COs9qr - 2x) , 
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Fig. 2. Proton-kaon correlation functions are shown for a 4 x 4 x 8 fm Gaussian 
source. Quantum effects mute the sensitivity predicted from classical Coulomb 
trajectories. This especially true at low Q = 2q. In the right-side panel, the 
correlation functions are scaled by to counter the l/Q^ fall-off. Incorporating 
strong interactions into the quantum calculations also affects the correlation. 



some of the sensitivity for qr/?i ~ 1 This especially harmful for small reduced 
masses, such as tttt correlations, since the correlations are smaller and need to be 
analyzed at smaller q. In fact, for cases where the Bohr radius is much larger than 
the source size (oq = l/z^e^ = 390 fm for tttt), the effect of Coulomb is a simple 
Gamow factor which is independent of r and thus provides no information. Figure 
Elshows correlations pK~^ correlations for a 4 x 4 x 8 fm Gaussian source. Compared 
to classical correlations quantum correlations are somewhat muted, especially the 
shape information which is determined by looking at the dependence on the direction 
of q. The muting is small for larger q, but there the correlation is small. Nonetheless, 
correlations are strong enough near 100 MeV/c to determine shape characteristics 
of the source. The analogous calculations for Fig. |21can be performed for any pair. 
Any interaction can be exploited to determine shape. For instance pA interactions, 
which have no Coulomb, also provide sensitivity to shape through shadowing. 

In Figure 12 the anisotropy of the source was seen by plotting the correlation 
function as a function of cosflg, where 9q was measured relative to the elongated 
axis. For more complicated three-dimensional shapes, especially where the source 
may not have reflection symmetries, determining characteristics of the shape can be 
more difficult. Spherical and Cartesian harmonics can then be used to express 
the correlation function in terms of a few coefficients at small i. These expansion 
coefficients are themselves functions of q and can be related to expansion coefficients 
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for the source functions. 

Ri,m{q) ^ J 4:TTr'^dr JCi{q,r) Se^m{r), (3) 

where these quantities are related to the correlation function, wave function and 
source distribution by the following definitions, 

= (4^)-V2 f dQgRiqjY.^U^g), (4) 



)Ce{q,r) = ^ J dcosOqr [\(j){q,r,cos9qr)f - l] Pi{cos9qr). 

The convenience of this expression is that the expansion coefficients labeled £ and 
m are related on a one-to-one basis with the corresponding coefficients in the source 
function. 

The disadvantage with using spherical harmonics is that it can be difficult to 
physically visualize the distortion associated with a given i and m, particularly for 
£ > 2. This can be remedied by using Cartesian harmonics. Cartesian harmonics 
of order £ = £x + £y + £z are linear combinations of spherical harmonics of the same 
order. Each Cartesian harmonic is of the form, 

Af{n) — nl.^riy'nl' + lower orders of n^; , rt^ , . (5) 

Some examples are given in Table ^ The disadvantage of Cartesian harmonics 
is that they are not orthonormal within a given £ multiplet. But, the {2£ + 1) 
harmonics with £x = 0,1 can be used to easily generate the other functions with 
the identity, Ae^ i„ i, = —Ae,.~2 l,+2 i, — Ai^-2 i„ e,-2- Since Cartesian harmonics 
can be expressed as a linear combination of spherical harmonics of the same £, 

R^q) = J Airr^dr lC,{q,r) 5^-{r), (6) 

where the coefficients are defined, 

R^q) ^ i^^±llll j ^Ric)AA^ql (7) 

S^q) ^ (^i±llll| ^5(q)^-(f7,). (8) 

References [El [71 provide many properties of Cartesian harmonics and gives expres- 
sions for orthonormality, for transforming to spherical harmonics, and for generating 
moments. 

Results for various moments were generated for p7r correlations from a blast- 
wave model and are displayed in Fig. |3| The blast-wave temperature was chosen to 
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Table 1. Cartesian harmonics for £ <= 4. Other harmonics can be found by 
swapping indices on both sides of the equation, e.g., a; <-> y. For example, given 
-^2,1,0 = nlriy - ny/5, swapping y ^ z gives A2,oa = nluz - 71^/5. 



■^1,0.0 — nx 


-4i,i,i 




^2,0,0 - nl - (1/3) 


^^4,0,0 = n^ — 


i6/7)nl + (3/35) 


A, 1,0 = n^ny 


^3,1,0 = nl 


Uy - i3/7)nxny 
7)nl-{l/7)nl + {l/35) 


^3,0,0 = ^1- (3/5)nj; 


.42,2,0 = nlnl - (1/ 
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Fig. 3. Decompositions in Cartesian harmonics of p7r"^ source and correlation 
functions are shown for a blast- wave source (blue circles). A Gaussian source (Red), 
chosen to reproduce the integrated moments, does a poor job reproducing the source 
and correlation functions for some moments. 
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be 120 MeV, the radius was 13 fm and the coUective velocity was 0.7c. Since hghter 
particles tend to originate from deeper inside the blast wave, the source functions 
have strong dipole moments, which are manifest in the £x — ^ Cartesian projection, 
where x refers to the outward direction. The £ = 2 moments describe the elliptic 
anisotropics, whereas the i — 1 terms can be more directly identified with the offset 
of the center of the ellipse. The {£x = l,£y = 0,iz = 2) term is driven both by 
the combination of the offset and elliptic asymmetries and by the boomerang shape 
caused by the inside-outside geometry of the cascade. The non-Gaussian nature of 
the shape is seen by comparison with the best-fit Gaussian source. Gaussians are es- 
pecially bad for fitting the behavior of the source function at larger r. Furthermore, 
the one-to-one correspondence for labels of the correlations and source function 
makes source imaging, which has been applied for angle-integrated correlations [ 
IHl El d|i tenable for higher harmonics. 

One of the great benefits of analyses based on spherical or Cartesian harmonics 
is that the shape can be studied as a function of r. For instance, at RHIC the 
source functions should have exponential tails along the beam axis due to the boost- 
invariant nature of the expansion, and perhaps in the outward direction due to 
resonances. The higher harmonics should then be pronounced at large r due to the 
decreasing justification of Gaussian fits. This is apparent in Fig. O 

Analyses such as those illustrated here can be performed for nearly every species 
measured at RHIC. For instance, with the last data run at RHIC, there should be 
enough data to analyze Ett correlations. There are dozens of possible analyses. In 
any model incorporating thermalization and collective fiow, source sizes and shapes 
for different species should follow simple systematic trends driven by the particle 
masses. Observing such trends would be of tremendous importance in solidifying 
our understanding of collision dynamics at RHIC. 
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